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TECHNICAL  NOTE  1925 


LINE -VORTEX  THEORY  FOR  CALCULATION  OF 
SUPERSONIC  DOWNWASH 

By  Harold  Mirels  and  Rudolph  C.  Haefeli 


SUMMARY 

The  perturbation  field  induced  by  a  line  vortex  in  a  super¬ 
sonic  stream  and  the  downwash  behind  a  supersonic  lifting  surface 
are  examined  for  the  purpose  of  establishing  approximate  methods 
for  determining  the  downwash  behind  supersonic  wings. 

Lifting-line  methods  are  presented  for  calculating  super¬ 
sonic  downwash.  An  unbent  lifting  line  (horseshoe -vortex  system) 
is  used  to  compute  the  downwash  behind  rectangular  and  triangular 
wings  and  the  results  are  compared  with  the  exact  linearized  solu¬ 
tions.  The  chordwise  position  of  the  lifting  line  giving  the 
best  average  agreement  with  the  exact  solution  is  noted  for  each 
wing.  A  bent  lifting  line  is  used  to  approximate  the  triangular 
wing,  and  the  results  are  in  good  agreement  with  the  exact  solu¬ 
tion  except  for  points  within  l/2  chord  of  the  wing  trailing  edge. 
The  use  of  a  bent  lifting  line  seems  promising  for  obtaining 
accurate  estimates  of  the  downwash  behind  swept  wings. 


INTRODUCTION 

Several  methods  for  obtaining  the  downwash  behind  supersonic 
wings  based  on  linearized  theory  have  been  presented.  These 
methods  utilize  conical  superposition  (reference  l),  doublet  dis¬ 
tributions  (references  2  and  3),  or  vortex  distributions  (refer¬ 
ences  4  and  5).  Each  of  these  methods  has  certain  disadvantages. 
Conical  superposition  is  restricted  to  wings  having  plan  forms 
composed  of  straight-line  segments  and  is  cumbersome  for  other 
than  trapezoidal  or  triangular  plan  forms.  The  doublet  and  vortex 
distributions  apply  to  arbitrary  plan  forms,  but  provide  integral 
expressions  for  downwash  that  are  generally  very  tedious  to  eval¬ 
uate  .  The  complexity  of  these  expressions  indicates  that  there  Is 
a  need  for  a  straight-forward  procedure  for  obtaining  reasonably 
accurate,  if  not  exact,  downwash  solutions. 
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A  logical  approach  to  the  development  of  a  simplified  super¬ 
sonic  downwash  theory  is  to  derive  the  supersonic  analogues  of  the 
line -vortex  procedures  that  have  proved  valuable  in  subsonic  theory. 
Certain  differences  exist,  however,  between  the  properties  of 
vortices  in  a  supersonic  stream  and  vortices  in  a  subsonic  stream. 
Similarly,  the  downwash  fields  behind  subsonic  and  supersonic  wings 
differ  in  certain  respects.  These  differences  must  be  investigated 
before  an  extension  of  subsonic  techniques  is  possible. 

The  present  report,  prepared  at  the  NACA  Lewis  laboratory,  has 
three  main  objectives:  (1)  The  downwash  field  induced  by  a  super¬ 
sonic  line  vortex  of  constant  slope  is  derived  and  discussed; 

(2)  the  downwash  behind  a  supersonic  lifting  surface  is  examined 
and  related  to  the  downwash  field  induced  by  line  vortices;  and 

(3)  lifting-line  methods  for  computing  downwash  are  presented  and 
calculations  based  on  these  methods  are  compared  with  the  exact 
linearized  solutions.  Zero-thickness  wings  (lifting  surfaces)  are 
considered  throughout. 


SYMBOLS 

The  following  symbols  are  used  throughout  this  report: 


X,Xi 


y,yt  >  Cartesian  coordinate  system 


z,z± 


perturbation  velocities 


J  /  unit  vectors 
k 


components  of  vortex  vector  (three-dimensional  field) 
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components  of  vortex  vector  (vortex  sheet) 


wing  span 

function  used  in  evaluation  of  finite  part  of  divergent 
integral  (equation  (7)) 

local  wing  pressure  coefficient 

effective  chord  of  bent  lifting  line 

root  chord  of  wing 


complete  elliptic  integral  of  second  kind  with  modulus 

/.  /b\2 


Gi(m)  = 


(Yt  -  m1X1)Q2m1Y1  -  Xj) 
ri  [<Yi  “  “i*!)2  +  (1  ‘  P2^2^2 


0±(0)  = 


XiYi 


ri  (Yi2  +  Zi2) 


spanwise  lift  distribution 
free -stream  Mach  number 
sD  ope  of  line  vortex  or  lifting  line 
^X±2  -  P2Yi2  -  P%2 

free -stream  velocity  (taken  in  x -direction) 


Xi  = 


x  -  x, 


y  -  y± 


z  -  z. 
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angle  of  attack 
cotangent  of  Mach  angle 
integration  interval 


>fip  -  1 


circulation 


perturbation-velocity  potential 

cpT  -cps 

free -stream  density- 
cross -sectional  area  of  vortex  tube 
a panwise  distribution  of  wing  circulation 
wing  circulation  at  midspan 

l/2 

resultant  vorticity  (|2  +  T|2  + 
finite  part  of  divergent  integral 
line  integral  along  closed  curve 


Subscripts : 

a,b  points  of  intersection  of  forward  Mach  cone  with  line  vor¬ 
tex,  lifting  line,  or  edge  of  vortex  sheet 

B  bottom  surface  of  z  =  0  plane 

C  line  integration 

i  =  1,2,3  .  .  .n  points  on  vortex  lines 

l  plan-form  leading  edge 

S  surface  integration 

T  top  surface  of  z  =  0  plane 

t  plan-form  trailing  edge 

Y  volume  integration 

o  variable  of  integration 
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Superscripts : 

-  value  of  function  at  point  of  discontinuity  when  approach¬ 

ing  from  negative  y-direction 

+  value  of  function  at  point  of  discontinuity  when  approach¬ 

ing  from  positive  y-direction 


THEORY 

General  Vortex-Field  Relations 

The  equations  relating  velocity  and  vorticity  distributions 
in  a  slightly  perturbed  supersonic  stream  are  derived  in  refer¬ 
ences  4  and  5.  These  relations  are  summarized  in  the  following 
section.  The  velocity  field  is  assumed  to  consist  of  a  major 
supersonic  free -stream  velocity  U  (taken  in  the  positive 
x-direction)  plus  small  perturbation  velocities  u,  v,  and  w 
such  that  the  linearized  equations  of  motion  are  applicable. 

Three-dimensional  vortex  field.  -  The  velocity  field  in  vector 
form  is 


q  =  (U  +  u)  i  +  vj  +  wk  (l) 

The  vortex  vector  field  to  is  defined  as  the  curl  of  the  velocity 
vector.  Thus 


(0  =  V  X  q 

=  |T  +  T)  j  +  £k  (2) 


where 
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Vortex  lines  are  lines  that  are  tangent  at  all  points  to  the  local 
vortex  vector  and  are  determined  from  the  equation 

*£  _  _  £z  /3\ 

I  -  q  t  (3) 

Reference  4  presents  the  following  integral  expressions  for 
the  perturbation  velocities  induced  by  the  three-dimensional  vortex 
field: 


u 


£ 

2k 


1 7  (Zq1  -  Yo» 

lJ V  r* 


dx  dy  dz^ 
o  ^  o  o 


(4a) 


where 


-m 


i2  r  (xo5  - 

v  =  -ir-l  I  — - 3-^—  dxD  dyD  dz0 


>2  r  (Y  I  "  X  *) 


r  r '  o*>  ~o'1'  ,  ,  , 

v  "  -  27  J  - 73 - dXo  *2, 


V  r° 


XQ  =  x  -  xD 


y0  -  y  -  y, 


zo  =  2  -  zo 


(4b) 


(4c) 


p2y  2 


2 

o 


The  subscript  o  indicates  a  variable  of  integration.  The  inte¬ 
gration  is  conducted  over  the  volume  V  included  in  the  forward 
Mach  cone  from  the  point  x,y, z .  That  is, 


j. 


i 
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The  symbol  /  designates  the  finite  part  of  a  divergent  inte¬ 

gral.  (See  references  4,  6,  or  7.)  The  procedure  for  obtaining 
the  finite  part  is  systematized  in  reference  7  as  follows: 

The  integrals  in  equations  (4a)  to  (4c)  are  of  the  form 


I 


A(3Cp) 

(x2  -  x0)3/2 


(5) 


The  upper  limit  in  equation  (5)  corresponds  to  limits  on  the  Mach 
cone  in  equations  (4a)  to  (4c).  The  finite  part  of  this  integral 
is 


where 


A(xo) 


(*2  -  *0)V2  ^ 


C  =  lim 
xo — *z2 


2A(x2) 

Vx2  "  xo 


-  J(x0) 


(6) 


(7) 


The  term  j(xQ)  is  the  indefinite  integral  of  equation  (5)  and 
J(xi)  is  the  indefinite  integral  evaluated  at  the  lower  limit. 

The  justification  for  this  procedure  is  presented  in  references  4, 
6,  and  7.  It  can  he  shown  that  finite  parts  are  treated  in  a 
manner  similar  to  that  for  ordinary  integrals.  The  rales  of  addi¬ 
tion,  differentiation  under  the  integral  sign,  transformation  of 
variables,  and  integration  by  parts  apply. 

Vortex  sheet.  -  If  vorticity  exists  as  a  surface  of  velocity 
discontinuity  in  the  zD  =  0  plane,  then  outside  this  plane,  |, 
rj,  and  £  are  all  zero,  but  in  the  zD  =  0  plane  £  =  0  whereas 
|  and  q  are  infinite.  However,  the  limits 


Y 
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=  lim  £  dzD 

dz0— ^0 

rj'  =  lim  T|  dzc 

dz0 — >0 

are  finite  and  are  given  by  (reference  5) 


The  subscripts  T  and  B  designate  velocities  on  the  top  and 
bottom  surfaces  of  the  z0  =  0  plane,  respectively.  The  perturba 
tion  velocities  induced  by  the  vortex  sheet  are  obtained  by  sub¬ 
stituting  equation  (8)  in  equations  (4a)  to  (4c).  In  particular, 
the  vertical -perturbation-velocity  field  (upwash)  is  given  by 


✓ 


w 


£ 

2jt 


To  5’  -V>’ 


tar,  d-y, 


(9) 


The  area  of  integration  S  includes  all  the  vorticity  in  the 
forward  Mach  cone  from  x,y,z. 

line  vortex.  -  The  vortex  lines  through  all  points  on  an 
infinitely  small  closed  curve  bound  a  vortex  "tube."  The  circu 
lation 


K  =  (OCT 


l/2 

(where  (0  =  (|2  +  T)2  +  £2)  is  the  resultant  vorticity  and  CT 
is  the  infinitesimal  cross-sectional  area  of  the  tube)  is  constant 
at  all  points  along  the  tube.  The  vortex-vector  components  at 
any  point,  in  terms  of  the  differential  distance  d lQ  along  the 
vortex  tube,  are 


£  = 


(0 


taD 

dl0 


v 


f 
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The  elemental  volume  is  dx0dy0dz0  =  GdZ0.  A  line  vortex  is  gen¬ 
erated  by  allowing  the  cross  section  of  the  vortex  tube  to  approach 
zero  while  maintaining  too  constant.  The  upwash  induced  by  such 
a  line  vortex  is  found  by  substituting  the  preceding  expressions 
into  equation  (4c)  and  equals 


K(YC  dx0  -  X0  dyQ) 


(10) 


where  the  integration  is  conducted  along  the  portion  of  the  line 
vortex  within  the  forecone  from  x,y,z.  A  line  vortex  cannot 
terminate  within  a  fluid  flow  field  but  must  either  form  a  closed 
curve  or  extend  to  infinity  or  a  boundary  of  the  field. 


Upwash  Induced  by  Line  Vortices 

Complicated  velocity  fields  can  be  generated  by  the  linear 
superposition  of  relatively  simple  fields.  It  will  therefore  prove 
useful  for  subsequent  developments  to  determine  the  upwash  field 
induced  by  line  vortices  of  constant  slope. 

Line  vortex  of  constant  slope  and  strength.  -  The  upwash  at 
x,y,z  due  to  a  line  vortex  of  constant  slope  m  and  strength 
K,  intersecting  the  forward  Mach  cone  at  x^y^  (fig.  1)  is 
(from  equation  (10)) 


_ _ _ m2  (y-mx)dy0 _ 

[(1-P2m2)y02  +  2m(P2my-x)y0  +  m2(x2-p2y2-P2z2)]  3/^2 

(11) 


The  integration  is  performed  in  appendix  A  and  yields 


£ 


w  = 


2  rt 


(y  -  xnx)  (ft2my  -  x) 


(s2  -  PV 


P2Z2) 


1/2 


[(y  -  mx)2  +  (1  -  P2m2)z2] 


(12) 
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The  finite  part  of  the  integral  in  equation  (11)  is  obtained  by 
substituting  the  lower  limit  into  the  indefinite  integral  (appen¬ 
dix  A).  No  contribution  appears  from  the  upper  limit.  By  a  trans¬ 
lation  of  coordinates,  the  upwash  due  to  a  line  vortex  from  iq^yq 
that  intersects  the  forward  Mach  cone  at  Xb,yt>  (fig*  2(a))  is 
found  to  be 


w 


(Ya  -  mXi)(p2nffi  -  Xj) 
rx  [(Yq  -  mXq)2  +  (1  -  P2m2)  z2 ' 


(13) 


where 


Xi  *=  x  -  xi 
=  7  -  7l 


rx  ^Xq2  -  PV  -  P2z2 

By  superposition  (fig.  2)  of  a  positive  line  vortex  K  from 
xq,yq  and  a  negative  vortex  -K  from  *2>72>  the  upwash  due  to 
a  line -vortex  segment  not  intersecting  the  forward  Mach  cone  is 


w  =  &T  [°2(m)  ' 


(14) 


where  the  notation 


Gi(m) 


ri 


(Yj  -  m1X1)(P2m1Y1  -  Xg) 

(Yi  -  mgXq)2  +  (1  -  P2mj2)  Zg2' 


(15) 


is  used.  The  subscript  for  G  indicates  the  appropriate  sub¬ 
scripts  for  X,  Y,  Z,  r,  and  m. 

Equation  (14)  applies  for  any  line  vortex  of  constant  strength 
and  slope.  The  circulation  K  is  positive  when  the  vortex  vector 
is  in  the  direction  of  integration.  When  the  line  vortex  inter¬ 
sects  the  forecone,  the  corresponding  limit  (infinite  Gg(m)  term) 
is  neglected  by  application  of  the  finite -part  concept. 


line  vortex  of  constant  slope  but  varying  strength.  -  Line 
vortices  may  coincide  to  form  a  resultant  line  vortex  whose  cir¬ 
culation  is  the  sum  of  the  strengths  of  the  superposed  components. 


» 


> 


i 
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The  general  line  vortex  is  then  one  of  variable  strength  K  along 
the  line  xQ  =  x0(y0)  having  the  local  slope  m^  =  dy0/dx0.  The 
upwash  induced  by  the  segment  from  X]_,y^  to  xg,y2  “ay  t®  written 


K(Y0  -  m^o) 


<iy. 


m  r  ' 
o  o 


(16) 


The  evaluation  of  equation  (16)  is  generally  tedious.  For  the 
particular  case  of  constant  slope,  however,  this  equation  may  be 
integrated  by  parts  to  yield  a  useful  expression  for  upwash.  Inas¬ 
much  as 


.  P2(Yq  -  “qXq)  _*M*) 
Vo3  d7° 

the  integration  by  parts  gives 


The  term  K0  [G0(m)]  is  disregarded  at  a  limit  corresponding  to  a 
point  on  the  forward  Mach  cone. 

Equation  (17)  is  a  generalization  of  equation  (14)  to  account 
for  variations  in  strength  along  a  line  vortex  of  constant  slope. 

Both  equations  are  of  fundamental  importance  because  appropriate 
distributions  of  such  line  vortices  will  be  used  to  simulate  a 
supersonic  lifting  surface. 

Characteristics  of  upwash  field  due  to  supersonic  line  vortex.  - 

1.  Infinite  line  vortices  of  constant  strength  and  slope.  An 
infinite  line  vortex  inclined  supersonically  ( | pm|  >1 )  to  the  free 
stream  is  shown  in  figure  3(a) .  Because  both  limits  in  equation  (14) 
are  neglected,  the  upwash  induced  by  this  line  vortex  is  zero.  This 
result  agrees  with  the  indications  of  oblique -airfoil  theory  because 
the  perturbation  velocities  are  zero  downstream  of  the  envelope  of 
the  Mach  cones  from  the  trailing  edge  of  a  two-dimensional  airfoil 
inclined  supersonically  to  the  free  stream. 
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The  inf ini te -line  vortex  along  the  line  (yD  -  yi)  ■  m(x0  -  xj) 
(fig.  3(h)),  inclined  subsonically  ( | Pm|<l),  has  one  limit  that 
intersects  the  forward  Mach  cone  whereas  the  other  extends  to 
infinity.  The  upwash  is  then,  for  0^{3m<l, 


w 


lim  G0(m) 
x0— »-• 
y« — > -® 


K  tj  1  -  82m2  (Yx  -  mix) 

2jc  (Yx  -  mXi)2  +  (1  -  P2m2)  z2 


(18) 


where  xi,yi  is  a  point  on  the  line  vortex.  The  upwash  is  infinite 
along  the  line  vortex.  For  m  =  0,  equation  (18)  becomes 


w  = 


(19) 


which  is  identical  with  the  expression  for  upwash  due  to  an  infinite 
vortex  parallel  to  the  stream  (along  the  line  y  =  y^)  in  an  incom¬ 
pressible  field. 

These  results  indicate  that  the  behavior  of  the  infinite  super¬ 
sonic  line  vortex  for  |Pm|>l  is  completely  different  from  that  of 
the  incompressible  flow  vortex.  However,  when  |pm|<l,  both  vortices 
have  similar  upwash  properties  in  the  vicinity  of  the  vortex  line 
and  are,  in  fact,  identical  for  m  ■  0. 


2.  Bent  line  vortices.  The  upwash  due  to  a  bent  line  vortex 
(fig.  3(c))  of  constant  strength  is 


w  =  [plfa")  "  Gl(m+)]  (2°) 

where  m^“  and  designate  the  slopes  of  the  line  vortex  before 

and  after  the  bend  at  xi,y^.  The  term  r^  appears  as  a  factor  in 
the  denominator  of  equation  (20)  so  that  the  upwash  exists  only  in 
the  aftercone  from  xj_,y-j_  and  is  infinite  on  the  cone  surface 
(except  in  the  z  =  0  plane ) .  This  infinite  value  of  upwash  is 
not  to  be  confused  with  the  infinity  introduced  at  the  intersection 
of  the  line  vortex  with  the  forward  Mach  cone,  which  is  eliminated 
by  application  of  the  finite -part  concept. 
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Linearized  Supersonic -Wing -Theory  Relations 

The  perturbation  velocities  on  supersonic  lifting  surfaces 
(zero-thickness  wings)  have  been  evaluated  for  a  large  variety  of 
plan  forms.  (See,  for  example,  references  7  to  9.)  Those  results 
will  now  be  utilized  to  determine  the  vortex  field  generated  by  a 
lifting  surface. 


Velocity  potential.  -  If  the  boundary  conditions  for  a  lifting 
surface  are  specified  in  the  z  =  0  plane,  the  u  and  v  veloc¬ 
ities  are  antisymmetric  and  the  w  velocities  are  symmetric  about 
this  plane.  Thus  the  velocities  on  the  top  and  bottom  surface  of 
the  z  =  0  plane  are  related  by 


'B 


(21) 


The  discontinuities  in  u  and  v  constitute  a  vortex  sheet. 
Because  the  flow  is  everywhere  irrotational,  except  across  this 
sheet,  a  perturbation-velocity  potential  cp  can  be  so  defined  that 


Sep 


aep  =  dx  +  5-  dy  +  ^  dz 


=  udx  +  vdy  +  wdz 


(22) 


When  the  undisturbed  flow  field  upstream  of  the  wing  is  considered 
to  be  of  zero  potential,  and  the  boundary  condition  requiring  that 
Urj>  =  u-g  =  0  off  the  wing  (antisymmetry  of  u  and  zero  lift  off 
wing)  is  applied,  the  potential  in  the  z  =  0  plane  may  be  obtained 
by  integrating  along  lines  of  constant  y. 


<PT 


*Pb 


(23) 


x 
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where  xj  is  the  equation  of  the  leading  edge  as  a  function  of  y. 
From  equation  (23)  it  may  be  concluded  that 

1.  Everywhere  in  z  ■  0  plane  except  behind  the  wing  leading 
edge  cpT  *  cpB  »  0. 

2.  At  a  given  span  station  9r  and  9b  remain  constant  for 
all  values  of  x  downstream  of  the  trailing  edge. 

Lilies  of  constant  potential  for  rectangular  and  triangular  wings 
are  shown  in  figure  4. 

Vortex  lines.  -  The  equation  for  the  vortex  lines  (from  equa¬ 
tion  (3))  is 


T) 1  dx  -  I'  dy  =  0  (24) 

When  the  following  expressions  (equations  (8),  (21),  and  (22))  are 
substituted 


V  -  -<TT  "  VB>  “  "2  “  2  3iT 


n  - 


ub  =  2  3T  =  "2  3T 


the  equation  for  the  vortex  lines  becomes 


8%,  8% 

^rdx  +  ^~dy  =  ^-dx  +  ^-dy  =  ° 


(25) 


Comparison  with  equation  (22)  shows  that  equation  (25)  represents 
lines  of  constant  potential.  Thus  the  vortex  lines  coincide  with 
the  lines  of  constant  potential  in  the  z  =  0  plane. 


Circulation.  -  The  circulation  included  between  two  points 
x1,y1  and  X£,y2  on  a  wing  is  given  by 


K 


=  ®udx  +  vdy  +  wdz 


(26) 


The  path  of  integration  is  arbitrary  except  that  the  path  should 
cross  the  z  =>  0  plane  only  at  the  two  specified  points.  If  the 
integral  is  taken  along  the  top  and  bottom  surfaces  of  the  z  =  0 
plane. 
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*2*2  *1>71 

K  =  /  (uq,  cbc  +  vrj.  dy)  +  I  (uB  dx  +  vB  dy) 

xl»yi  ^x2^2 


-  <**,2  -%,!>  +  (<%,!  -f>B,2)  *  *%  -  *Pl 


(27) 


where  Acp  equals  cp^  -  cpjj  and  represents  the  jump  in  potential  at 
the  point.  The  quantity  Acp  is,  in  fact,  the  doublet  potential 
(reference  2).,  so  that  the  net  circulation  between  two  points  equals 
the  difference  in  the  doublet  potential  between  those  points.  The 
equivalence  of  a  doublet  distribution  and  a  vortex  distribution 
indicates  that  the  flow  about  a  lifting  surface  can  be  calculated  cn 
either  basis. 


Circulation  ana  lift.  -  The  lift  per  unit  span  is  given  by  the 
chordwise  integration 


B(y) 


(c 


P>B 


“  Cp>T)  dx 


(28) 


After  substitution  of  the  linearized  values  for  pressure  coefficient 

2  Urn  2 

CP,T 

2uB  g 

°p,B  =  -  — - U  35T 

and  integration,  equation  (28)  becomes,  because  Aq>j  =  0, 

L(y)  =  pU  A (29) 

The  factor  Atp^  is  the  circulation  included  between  the  leading 
and  trailing  edges  at  the  spanwise  station  tinder  consideration. 

When  this  circulation  is  designated  T,  equation  (29)  becomes 


r 


L(y)  =  pur 


(30) 
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which  is  the  familiar  incompressible -flow  relation.  Also 

dy  dy 

=  (VT  “  TB^ 

=  (51) 

Equation  (31)  relates  the  shed  vorticity  to  the  rate  of  ohange  of 
spanwise  loading. 


APPLICATIONS  TO  CALCULATIONS  OF 
SUPERSONIC  DOWNWASH 

The  vertical  perturbation  velocities  due  to  a  supersonic  line 
vortex  of  constant  slope  have  been  presented  in  a  form  that  permits 
analytical  or  mechanical  evaluation  (equations  (14)  and  (17)).  The 
vortex  distribution  associated  with  a  lifting  surface  has  also  been 
discussed.  These  relations  will  be  used  to  develop  exact  and  approx¬ 
imate  methods  for  calculating  downwash  behind  lifting  surfaces  at 
supersonic  speeds. 


Downwash  an  Infinite  Distance  behind  Wing 

The  vertical -perturbation -velocity  field  behind  a  supersonic 
lifting  surface  (from  equation  (9))  is 


w 


FT  do*' -v i'> 

/4 


dxD 


dyD  - 


*ot'  ,  . 

—  3-  ^0  d7o 
ro 


(32) 


where  the  integration  over  the  plan  form  is  designated  by  Sp  and 
over  the  wake  by  Sy.  As  x  becomes  infinite,  X0  also  becomes 
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infinite  in  the  integral  for  the  hound  (plan  form)  vortices.  This 
integral  then  "becomes  zero  because  Xq  is  of  higher  order  in  the 
denominator  than  in  the  numerator.  Thus,  only  the  integration  over 
the  trailing  vortex  sheet  contributes  to  the  vertical  perturbation 
velocities  at  infinity.  The  trailing  vortex  sheet  may  be  con¬ 
sidered  to  consist  of  elemental  vortices  of  infinite  length  along 
y  =  constant  lines,  each  having  the  strength 


4K  =  I'  cLy0  = 

»,y,z  due  to 
tion  (19) 


ar 


dy0.  The  vertical  perturbation  velocity  at 


the  elemental  vortex  along  y  *  yQ  is,  from  equa- 


so  that  the  vertical -perturbation-velocity  field  at  infinity  is 
given  by 


(33) 


where  b/2  is  the  semispan.  The  velocity  field  obtained  from 
equation  (33)  is  identical  to  that  induced  by  a  subsonic  wing  with 
the  given  span  loading.  The  velocity  field  at  infinity  is  thus 
independent  of  Mach  number  (excluding  the  influence  of  Mach  number 
on  dr/dy0) .  This  result  has  been  derived  in  references  1  and  2 
by  other  methods.  The  evaluation  of  equation  (33)  is  relatively 
simple  and  may  be  used  to  approximate  the  downwash  several  chords 
behind  a  supersonic  wing. 


Regions  of  Infinite  Downwash 

Approximate  solutions  may  modify  or  introduce  singularities 
in  the  downwash  field.  It  is  therefore  of  interest  to  establish 
the  regions  for  which  linearized  theory  indicates  an  Infinite 
downwash. 

Infinite  downwash  in  z  »  0  plane  ♦  -  The  vertical  perturba- 
batlon  velocities  in  the  z  =  0  plane  that  exist  an  infinite 
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distance  "behind  a  wing  having  a  discontinuity  in  dr/dy0  will 
first  he  considered.  These  velocities  are  determined  from  the 
equation 


b 


2 


(34) 


The  Cauchy  principal  value  is  required  for  points  on  the  vortex 
sheet.  For  the  particular  case  of  triangular  loading  (fig.  5(a)) 


» 


dr  _ 
4y0  =  b 


in  the  integration  interval 


T-yo-°  and 


dr  _  srm 

iy0  =  ”  * 


in  the  interval  0<  yQ<  j 

dp 

When  these  values  for  j— 

dyc 

integral  yields 


(rm  is  the  circulation  at  the  midspan) . 
are  substituted  in  equation  (34),  the 


w  = 


(35) 


Infinite  upwash  exists  along  the  lines  y  =  dtb/2  whereas  infinite 
dcrwnwash  exists  along  the  line  y  *  0.  These  infinite  values  are 
due  to  the  discontinuities  in  the  spanwise  vorticity  distribution 
and  apply  for  all  points  on  these  lines  downstream  of  the  wing 
trailing  edge. 

In  general,  if  (dT/dy0)"  and  (dT/dy0)+  represent  a  dis¬ 
continuity  in  the  rate  of  change  of  spanwise  loading  at  station 
y  =  y^,  infinite  vertical  perturbation  velocities  will  exist  along 
y  =  yj  downstream  of  the  trailing  edge.  For  (dr/dy0)“ < (dT/dy0)* 
infinite  upwash  will  exist,  and  for  (dr/dyQ)~  >(dr/dyQ)+  infinite 


* 
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downwash  will  exist  along  this  line.  Such  discontinuities  in 
dr/dyo  originate  "both  at  a  wing  tip  and  at  those  points  along  a 
supersonic  trailing  edge  where  the  plan-form  slope  is  discontinuous 
and  the  local -wing-pressure  coefficient  is  not  equal  to  zero.  This 
discontinuity  in  vorticity  may  he  verified  hy  the  methods  employed 
in  reference  1  for  finding  the  upwash  and  sidewash  directly  behind 
a  supersonic  trailing  edge. 

The  discontinuity  in  shed  vorticity  at  the  tips  of  an  ellip- 
tically  loaded  wing  is  a  special  case  of  the  previously  stated  rule 
governing  infinite  vertical  perturbation  velocities  in  the  z  »  0 
plane .  For  wing  loading  given  by 


r 


m 


the  shed  vorticity  is 


dr 

dy0 


y0 


-;jyo 


and  the  vertical -perturbation-velocity  field  at  infinity,  in  the 
z  =  0  plane,  is 


Integration  yields 


w 


for 
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and 


As  Indicated  in  figure  5(b),  the  vertical  perturbation  velocity  is 
discontinuous  at  y  *  ±b/2,  but  is  bounded  for  all  points  on  the 
vortex  Bheet. 

Infinite  downwash  on  Mach  cone  from  wing  tip.  -  Reference  3 
indicates  infinite  vertical  perturbation  velocities  on  the  down¬ 
stream  Mach  cones  from  the  tips  of  a  triangular  wing.  This  result 
will  now  be  extended  to  apply  to  any  wing  tip  formed  by  the  inter¬ 
section  of  a  subsonic  leading  edge  and  a  supersonic  trailing  edge 
provided  that  the  slope  of  the  subsonic  edge  is  not  zero  at  the 
tip  (fig.  6) . 

The  contribution  to  the  vertical -perturbation-velocity  field 
due  to  the  bending  of  an  elemental  vortex  at  the  trailing  edge 
(fig.  6)  is,  from  equation  (20), 

iy"  dFo  r  -j 

dw  =  -  ■■°-2T--  [_G0(m-)  -  Go(0)J  (36) 


where 


G„(0)  =  ±  _ 
0  <To 


*0*0 


+  ZoZ) 


The  vertical  perturbation  velocity  at  a  point  on  the  Mach  cone  from 
this  tip  due  to  the  bending  of  the  vortices  is  found  by  integrating 
(along  the  trailing  edge) 


w 


G0(nf)  -  Go(0) 


dr 


dy0 


<*y0 


(37) 
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Equation  (37)  In  the  expanded  form  becomes 


w  * 


(Yp  -  m0-X0)(P2m0-ro  -  Ip)  XQY0 

(Y0  -  Vxo)2  +  (1  -  P2^2)*8  +  Y02  +  z2 


1  dr  , 

ro  4y0  4y< 


(38) 


The  limits  of  Integration  are  roots  of  r0,  so  that  the  factor 
l/rQ  Is  singular  of  order  l/2  at  the  limits  aid  the  integral  is 
improper  (assuming  (m0”)yo  m  q  £  0).  The  convergence  of  the  integral 

depends  on  the  nature  of  dr/dyQ  at  the  limits.  Appendix  B  shows, 
however,  that  dT/dy0  is  also  singular  of  order  l/2  at  the  tip. 

The  combined  singularity  causes  equation  (38)  to  diverge  at  the 
upper  limit  and  results  in  logarithmically  infinite  vertical  per¬ 
turbation  velocities  on  the  Mach  cone  from  the  tip.  The  divergence 
is  a  consequence  of  both  the  singularity  in  dl/dy0  and  the  singu¬ 
larity  on  the  Mach  cone  from  a  bend  in  the  elemental  line  vortex. 


The  infinite  vertical  perturbation  velocities  on  the  Mach  cone 
from  the  tip  do  not  appear  in  the  z  =  0  plane  inasmuch  as  equa¬ 
tion  (38)  then  reduces  to 


v 


<_  V*o2  -  eV  ar  a 

(Y0  -  mo'Xo)Y0  dy0  ay° 


and  the  singularity  due  to  r0  is  no  longer  present. 


Approximate  Downvash  Solutions 

Several  approximate  methods  for  obtaining  downwash  were  con¬ 
sidered.  Methods  based  on  a  lifting  line  seem  the  most  promising 
and  are  discussed  in  the  following  sections .  A  bent  lifting  line 
is  proposed  for  determining  the  downwash  behind  a  swept  wing  and 
an  unbent  lifting  line  (horseshoe  vortex  system)  for  determining 
the  downwash  behind  an  unswept  wing.  These  methods  are  applied  to 
compute  the  downwash  behind  triangular  and  rectangular  wings  and 
the  results  are  compared  with  the  exact  linearized  solutions . 
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A  lifting  line  concentrates  the  chordwise  loading  into  a  line. 
Thus  the  bound  circulation  is  represented  by  a  line  vortex  of  vari¬ 
able  strength  (k  =  T  *  ,  whereas  the  trailing -vortex  system 

maintains  the  same  vorticity  (£'  =  -dT/dy0)  but  now  originates  at 
the  line  rather  than  at  the  trailing  edge. 


Bent  lifting  line.  -  A  lifting  line  approximating  the  section 
centers  of  pressure 


(°p,B  “  Cp>T)(x  -  xj)  dx 


xcp  ~  + 


/ 


(°p,B  “  Cp,T^  dx 


seems  to  be  a  reasonable  representation  of  a  sweptback  or  swept- 
forward  wing.  In  order  to  facilitate  dcwnwash  calculations,  the 
line  of  section  centers  of  pressure  can  be  approximated  by  two 
straight-line  segments,  each  connecting  an  end  point  to  the  mid¬ 
point  of  the  line.  The  result  is  a  bent  lifting  line  (fig.  7)  of 
span  b  and  effective  chord  ce  (x-distance  between  midpoint  and 
end  point b  of  lifting  line).  For  a  sweptback  lifting  line,  the 
slopes  at  the  midpoint  are 

(“o’Vo  ■  0  c  "  2c^ 

and 

(“o+V0  =  0  =  2c^ 

The  vertical  perturbation  velocities  induced  by  the  bound 
vortices  can  be  determined  from  equation  (17)  and  equal 


9 


« 


* 


where  the  integration  is  conducted  along  the  lifting  line  (figs.  8(a) 
and  8(b)).  The  vertical  perturbation  velocities  induced  by  the 
trailing -vortex  system  are  obtained  from 
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(40) 


The  term  Go(0)  is  as  defined  for  equation  (36)  and  the  integration 
is  again  conducted  along  the  lifting  line.  The  vertical  perturba¬ 
tion  velocities  induced  by  the  complete  lifting-line  system  (bound 
and  trailing  vortices)  are  then  the  sum  of  equations  (39)  and  (40) 
and  equal 


v  »  §*  [Go(m")  -  Go("i+)] 


y0 


=  0 


1_ 

2n 


^b 


G0(m) 


0o(0)] 


dT 

iyr 


4y, 


ya 


(41) 


The  value  nu  =  (m  ")  n  is  used  for  the  integration  interval 

°  y0  =  u 

ya<yo<0  and  =  (m0+)T  _  q  is  used  for  the  interval 

o^y0^yb*  ° 


Equation  (41)  is  the  expression  for  the  vertical  perturbation 
velocities  behind  a  swept  wing  using  the  bent -lifting-line  approxi¬ 
mation.  This  equation  can  also  be  considered  as  derived  from  a 
superposition  of  a  system  of  bent  line  vortices  of  constant  strength 
(equation  (20)),  as  indicated  in  figure  9.  The  nonintegral  term  of 
equation  (41)  is  then  the  contribution  of  the  bent  line  vortex  of 
strength  rm  whereas  the  integral  term  represents  the  contribution 

of  the  elemental  line  vortices  of  strength  S —  dyD . 

dy0  • 


The  integral  term  of  equation  (41)  may  be  evaluated  by  analyt¬ 
ical  or  mechanical  methods .  When  mechanical  methods  are  used,  the 
singularities  in  the  integrand  must  be  isolated.  Suitable  proce¬ 
dures  for  isolating  the  commonly  encountered  singularities  are  as 
follows : 

(1)  Singularity  due  to  intersection  of  forward  Mach  cone  with 
lifting  line.  The  integrand  in  equation  (41)  is  infinite  at  the 
intersection  of  the  forward  Mach  cone  and  the  lifting  line  (for 
z  /  0) .  The  contribution  6w  of  the  integral,  for  the  interval 
y-b  -  e-tj^yo^y-b  (fig.  8(a)),  to  the  vertical  perturbation  velocity 
field  may  be  written 


*  =  -  £r  LG<»  -  Go«»J  ro  g-  -  LGo(“)  -  Go(°)J  r0  g 
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for 


where  mQ  Is  the  slope  of  the  lifting  line  at  y-5.  A  similar  pro¬ 
cedure  applies  for  the  singularity  at  ya. 

(2)  Singularity  in  dr/dy0  at  wing  tip.  The  vertical  pertur¬ 
bation  velocity  at  a  point  whose  forward  Mach  cone  intersects  the 
edges  of  the  trailing  vortex  sheet  (fig.  8(b))  is  obtained  from 

equation  (41)  by  integration  between  the  limits  -  |i<y0<g>.  If 

dT/dy0  is  singular  at  these  limits,  the  singularity  may  be  isolated 
by  a  procedure  similar  to  that  used  in  equation  (42).  Thus,  the 
contribution  to  the  vertical  perturbation  velocity  field,  from  the 

■jj  *jj 

integration  interval  77  ~  eb  -  ^o-  "written 
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(3)  Singularity  at  points  on  vortex  sheet.  The  Cauchy  prin¬ 
cipal  value  of  the  integral  (equation  (41))  is  required  for  deter¬ 
mining  the  vertical  perturbation  velocity  at  a  point  on  the  vortex 

v  v 

sheet  z  =  0,-  |'<y<|-.  For  this  case,  considering  the  interval 

y  -  e  <  y  <y  +  e  (fig.  8(c)),  the  contribution  of  the  integral  to 
the  vertical -perturbation-velocity  field  is 


8w 


-  eV  dF 

(Y0  ”  mc?o^o 


(44) 


If  dr/dy0  can  be  approximated  by  the  first  two  terms  of  the 
Taylor's  expansion 
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and  e  is  sufficiently  small  that 


moVxo2  -  P2Yq2  s  -  1 

(Y0  ~  moX0) 

equation  (44)  may  he  written 


Inasmuch  as  dT/dy0  is  an  odd  function  in  y0  for  a  wing  symmet¬ 
rical  about  the  midspan,  dT/dy0  is  either  discontinuous  or  zero 
for  yD  =  0.  If  dr/dy0  is  discontinuous,  the  vertical  perturba¬ 
tion  velocities  are  infinite  along  this  line. 


Unbent  lifting  line.  -  The  unbent  lifting  line  (horseshoe  - 
vortex  system)  appears  to  be  a  reasonable  representation  for  an 
upswept  wing.  The  use  of  such  a  lifting  line  considerably  sim¬ 
plifies  the  calculations. 


The  vertical  perturbation  velocities  induced  by  an  unbent 
lifting  line  (along  the  y-axis)  are  obtained  from  equation  (41) 
by  setting  hIq  =  «  and  are  equal  to 


w 


x  Y0(rp2  -  e2z2)  ar 
r0(x2  -  P2z2)(Y02  +  z2)  *y0  70 


The  singularities  are  isolated  as  follows: 


(46) 
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(1)  Singularity  due  to  intersection  of  forward  Mach  cone  with 
lifting  line. 


5w  =  - 


yb-eb 


*  Y0(r02  -  eV)  dr  dr  \ 

(x2  -  82z2)(Y02  +  z2)  dyo  \x  Yo  &o)  ro 

0  y0=yb 


_i_  /p2z2  dr 

2rt  lx  Y0  dyD 


y0=yb 


yb-£b 


where 


2  «1n-i  gt 

/  "  F  1  v^b^ 

Jyb-€b 

(2)  Singularity  in  dF/dy,..  at  the  tip. 


x  YQ(ro2-p2z2) 


/b  |r0(x2-82z2)(Y02+z2)  rQ(x2-P2z2)(Y02+z2) 

2-S 


x  YQ(ro2-P2z2) 


b 

y0-2 


dr  ,,  l. 


1  1  ^0<r02  -  ^ 


2rt  r^x2  -  P2Z2)(YQ2  +  z2) 


(r)  b 

b  7o  “  2 


yQ  -  2 


(3)  Singularity  at  points  on  vortex  sheet 


c  /dzr 

“Wo* 


y0-y 
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The  chordwise  location  of  the  unbent  lifting  line  that  will  give 
the  best  average  agreement  with  the  exact  linearized  solution  is 
still  to  be  determined.  According  to  the  techniques  used  in  subsonic 
wing  theory,  an  unbent  lifting  line  at  the  wing  center  of  pressure 
should  be  a  good  first  approximation.  However,  further  investi¬ 
gation  is  required.  It  may  be  possible  to  determine  the  best  loca¬ 
tion  for  each  general  class  of  plan  forms  by  comparing  the  lifting¬ 
line  solution  with  the  exact  linearized  solution  for  several 
representative  plan  forms. 

Examples  and  discussion.  -  A  bent  lifting  line  and  an  unbent 
lifting  line  will  be  used  to  obtain  solutions  for  the  downwash  in 
the  y  =  0  plane  behind  triangular  wings  having  subsonic  leading 
edges.  The  chordwise  distribution  of  wing  loading  suggests  the  use 
of  the  bent  lifting  line.  The  unbent  lifting  line  will  be  used  for 
purposes  of  comparison.  An  unbent  lifting  line  also  will  be  used 
to  compute  the  downwash  in  the  y  =  0  plane  behind  rectangular 
wings.  The  chordwise  location  of  the  unbent  lifting  line  giving 
the  best  average  agreement  with  the  exact  linearized  solutions 
will  be  noted  for  both  the  triangular  and  the  rectangular  wings. 

In  order  to  simplify  the  expressions,  M  =  ^2  (that  is,  P  *  1) 
will  be  assumed. 


The  spanwise  circulation  distribution  for  a  triangular  wing 

of  span  b  and  root  chord  cr,  having  subsonic  ^~—<lj  lea&in8 
edges,  is  (from  equation  (24)  of  reference  2)  r 


r 


where  E  is  the  complete  elliptic  integral  of  the  second  kind 
with  modulus 

r  _  cfljb 
m  E 


-( 


_b 

2  c, 


The  circulation  at  the  midspan  is 


and  the  rate  of  change  of  circulation  is 


*■ 


* 


dr  4au  yc 
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The  expression  for  the  vertical  perturbation  velocities  in  the 
y  =  0  plane  using  the  bent-lifting-line  approximation  of  fig¬ 
ure  10(a)  is 


w 


[G0(m+)  -  Go(0)] 


d?o 


(47) 


where  (mr.+) 

0  yD=  0 


— .  Equation  (47)  is  obtained  from  equation  (41) 
cr 


and  the  properties  of  an  even  function.  The  vertical -perturbation- 
velocity  field  behind  an  unbent  lifting  line,  having  the  same 
loading  (fig.  10(b))  is  from  equation  (46), 


w  = 


x(r02  -  z2)  y02dy0 


/  2  2 w  2  2\ 

rc(x  -  z  )(y0  +  z  ) 


—  r2 
*2  7° 


(48) 


The  spanwise 
b 
c 


circulation  for  a  rectangular  wing  of  aspect 


ratio  — ?2  for  the  tip  region  ^  -  cr^y0f=;j-  is 


and  for 


(49a) 


r  =  2dU  cr  (49b) 

Equations  (49a)  and  (49b)  were  obtained  from  equation  (20)  of  refer¬ 
ence  8.  The  wing  is  illustrated  in  figure  10(c).  The  rate  of 

Id  Id 

change  of  spanwise  circulation  for  ^  -  cr<  yQ<  g-  is 
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The  vertical -perturbation-velocity  field  in  the  y  -  0  plane 
behind  an  unbent  lifting  line  having  this  loading  is 


w 


x  (rp2  -  z2)  yQ 
rQ(x2  -  z2)(yQ2  +  z2) 


(50) 


Equation  (47)  was  evaluated  to  determine  the  downwash  along 
the  line  y  =  0,  z  =  0  and  along  the  line  y  =  0,  z  =  0.1  b 
for  triangular  wings  with  aspect  ratios  (2b/cr)  of  1.6  and  3.2. 
Equation  (48)  was  evaluated  for  the  same  wings  at  points  along 
the  line  y  =  0,  z  =  0.  Equation  (50)  was  evaluated  at  points 
along  the  line  y  =  0,  z  =  0  for  rectangular  wings  with  aspect 
ratios  (b/cr)  of  2  and  4.  The  integrations  were  performed  mechan¬ 
ically  for  equations  (47)  and  (50)  and  analytically  for  equa¬ 
tion  (48).  The  results  are  compared  with  the  exact  linearized 
solutions  obtained  from  references  1  to  3  in  figures  11  to  14. 

The  bent -lifting -line  solutions  for  the  triangular  wings  are 
shown  in  figures  11  and  12.  The  discontinuity  in  the  curves  for 
the  z  =  0.1  b  plane  (figs.  11(b)  and  12(b))  designates  large 
negative  (upwash)  values  that  become  infinite  on  the  Mach  cone 
from  the  tip,  as  indicated  in  the  discussion  associated  with 
equation  (38) .  The  agreement  with  the  exact  solutions  is  good 
for  al 1  points  except  those  within  l/2  chord  from  the  trailing 
edge.  The  correlation  is  unexpected  because  of  the  large  contri¬ 
bution  of  the  nonintegral  term  of  equation  (47).  This  term  is 
associated  with  the  bending  of  the  lifting  line  at  the  midpoint. 
This  bend  was  artificially  introduced .  The  agreement  with  the 
exact  downwash  solutions,  however,  indicates  that  the  bent  lifting 
line  is  a  good  average  representation  of  the  triangular-type  wing 
(at  least  in  regard  to  the  downwash  solution  in  the  region  of 
the  line  y  =  0,  z  =  0) .  This  method  should  give  even  better 
agreement  when  used  to  represent  a  sweptbach  wing  such  as  the  one 
represented  in  figure  7,  because  the  bent  lifting  line  would  then 
more  closely  approximate  the  actual  vorticity  distribution. 
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The  uribe nt -1 i f t i  ng -1 ine  solution  for  the  triangular  wings  is 
presented  in  figure  13.  The  lifting  line  is  placed  at  the 
3/4-chord  point  to  give  the  best  average  agreement  with  the  exact 
linearized  solution.  (The  center  of  pressure  for  these  wings  is 
at  the  2/3-chord  point.)  The  agreement  is  not  as  good  as  that 
obtained  with  the  bent  lifting  line  and  indicates  that  the  bent 
line  is  more  suitable  for  computing  downwash  behind  triangular 
wings.  The  smaller-aspect -ratio  wing  (2b/cr  =  1.6)  is  in  better 
agreement  with  the  exact  linearized  solution  because  of  the  rapidity 
with  which  the  downwash  approaches  the  asymptotic  value  at  infinity. 

The  unbent -lifting -line  solution  for  the  rectangular  wings  is 
presented  in  figure  14.  The  downwash  obtained  with  the  lifting 
line  at  the  l/2 -chord  point  is  in  excellent  agreement  with  the 
exact  linearized  solution.  (The  location  of  the  center  of  pressure 
is  at  the  4/9-chord  and  10/21-chord  points  for  aspect  ratios  of  2 
and  4,  respectively.)  These  calculations  suggest  that  the  unbent 
lifting  line  is  well  suited  for  computing  the  downwash  behind 
unswept  supersonic  wings  and  that  the  best  chordwise  position  for 
the  unbent  lifting  line  is  at  or  slightly  downstream  of  the  wing 
center  of  pressure. 

A  more  accurate  estimate  of  the  downwash  in  the  z  =  0  plane 
close  to  the  trailing  edge  may  be  obtained  by  judiciously  fairing 
the  curve  obtained  by  the  line -vortex  method  to  the  known  value 
of  downwash  at  the  trailing  edge.  At  a  subsonic  trailing  edge 
satisfying  the  Kutta  condition,  -w/aU  is  unity;  whereas  at  a 
supersonic  trailing  edge,  -w/oflj  can  be  computed  by  the  method 
presented  in  reference  1. 

It  should  be  noted  that  linearized  theory  neglects  the  effect 
on  the  downwash  of  the  friction  wake  and  the  displacement  and  dis¬ 
tortion  of  the  trailing  vortex  sheet.  An  experimental  program, 
such  as  that  reported  in  reference  10  for  wings  in  subsonic  flight, 
is  ultimately  required  to  determine  the  necessary  modifications  to 
linearized  theory  that  will  result  in  good  agreement  between  theory 
and  practice. 


SUMMARY  OF  ANALYSIS  AND  APPLICATIONS 

The  perturbation  field  induced  by  a  line  vortex  in  a  super¬ 
sonic  stream  and  the  downwash  field  behind  a  supersonic  lifting 
surface  have  been  examined  for  the  purpose  of  establishing  approx¬ 
imate  methods  for  the  calculation  of  supersonic  downwash. 
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An  infinite  line  vortex  of  constant  strength  and  slope,  super¬ 
sonically  inclined  to  the  free  stream,  induces  no  perturbation 
field.  A  subsonically  inclined  line  vortex  has  properties  similar 
to  those  of  a  vortex  in  an  incompressible -flow  field.  Bends  in  a 
line  vortex  induce  infinite  vertical  perturbation  velocities  on 
the  surface  of  the  downstream  Mach  cone  from  the  bend  (except  in 
the  z  =  0  plane) . 

The  downwash  field  behind  a  supersonic  lifting  surface  differs 
from  that  behind  a  subsonic  wing  in  several  respects.  For  a  super¬ 
sonic  lifting  surface,  discontinuities  in  shed  vorticity  occur  at 
those  points  along  a  supersonic  trailing  edge  where  the  plan-form 
slope  is  discontinuous  and  the  local  pressure  coefficient  is  not 
zero.  These  discontinuities  lead  to  singularities  in  the  downwash 
field  in  the  z  =  0  plane.  Also,  the  vertical  perturbation  veloc¬ 
ities  are  logarithmically  infinite  on  the  downstream  Mach  cone 
from  a  wing  tip  formed  by  the  intersection  of  a  subsonic  leading 
plan-form  edge  and  a  supersonic  trailing  plan-form  edge. 

A  bent -lifting -line  method  has  been  proposed  for  the  solution 
of  the  downwash  field  behind  swept  wings.  When  applied  to  a 
triangular  wing,  this  method  gave  results  that  were  in  very  good 
agreement  with  the  exact  linearized  solution  for  points  near  the 
line  y  =  0,  z  =  0  except  for  points  within  l/2  chord  of  the 
wing  trailing  edge. 

An  unbent  lifting  line  (horse shoe -vortex  system)  has  been 
proposed  for  unswept  wings.  This  method  was  applied  to  determine 
the  downwash  behind  rectangular  wings  with  aspect  ratios  b/cr  of 
2  and  4.  Excellent  agreement  was  obtained  for  both  aspect  ratios 
by  placing  the  lifting  line  at  the  l/2 -chord  point. 


Lewis  Slight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  April  29,  1949. 
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APPENDIX  A 


UPWASH  INDUCED  BY  LINE  VORTEX 

The  upwash  Induced  "by  the  line  vortex  segment  of  figure  1  is, 

_ _ m2(y-mx)dy0 _ 

q  [(1  -P2n?)y02  +  2m(p2my-x)y0  +  n?(x2-P2y2-P2z2)J  3^2 

(11) 


The  upper  limit  is  at  the  intersection  of  the  line  vortex  (y0  =  mxQ) 
with  the  trace  of  the  forward  Mach  cone  in  the  z  =  0  plane 

[(x  -  xc)2  -  P2(y  -  yD)2  -  P2z2  =  o]  and  is  therefore  the  appropriate 
root  of  the  expression  appearing  in  the  denominator  of  the  integrand. 
The  roots  are 


-m( P2my  -  x)  ±  Pm  V(y  -  mx)2  +  (1  -  P2m2)z2 
=  (i  -  eV) 

so  that  equation  (11)  may  he  rewritten 


3 

fPr-jWi..  E. 


iy, 


p2m2K(y  -  mx)  U  £(ya  -  y0)(Xb  -  y0)] 

From  equations  (6)  and  (7), 


5/2 


fl  =  -  J(0)  -  C 


where 


J(0) 


2(ya  +  yb) 


(ya  “ 


(Al) 


(A2) 


(A3) 
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and 


C  =  lim 

y0 — ^yb 


=  lim 
Jo — ry^, 


=  0 


2A(yb)  w  \ 
T==  "  J(yo} 
p-b  -  Jo 


2(ya  +  J-o  -  2y0) 


, ,  J  , -  (ya-yt)2 

L(ya-yb)  f[j\>-J0 


Thus 


/I  -  - 


2(ya  +  y-b) 


(ya  "  y-b^V^b 


(A4) 


Substituting  equation  (A4)  in  equation  (A2)  and  solving  for  w 
yields 


w  =  - 


JL 

2« 


(y-mx)  ( p2my-x) 


X  _ 

/  2  „2  2  q2  2x2  f /  \2  /.  Q2  2\  2  1 

(x  -P  y  -P  z  )  [(y-mx)  +  (1-P  m  )z  J 


(A5) 
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APPENDIX  B 


LOADING  IN  VICINITY  OF  WING  TIP 

The  nature  of  the  loading  in  the  vicinity  of  a  wing  tip  formed 
"by  the  intersection  of  a  subsonic  leading  edge  with  a  supersonic 
trailing  edge  (fig.  6)  is  to  he  determined. 

According  to  linearized  theory,  uj»  in  the  region  of  a  sub¬ 
sonic  leading  edge  is  singular  of  the  form 


where  xj  is  the  equation  for  the  leading  edge  as  a  function  of 
y0.  This  relation  may  be  deduced  from  equation  (11)  of  refer¬ 
ence  10.  The  corresponding  wing  circulation  is 

f1* 

/  dxQ  _ _ 

r  =  Acp^.  a  2  /  "  ■-  =  4  a/X|.  -  Xj  (Bl) 

/  4xo~xl 
xl 


The  derivative  of  equation  (Bl)  is 


„  ^(xt  -  Xj) 

dr  dy0 

fa  “  xl 


(B2) 


where  (x^  -  xj)  is  the  wing  chord  as  a  function  of  yc.  However, 
(x^  -  xj)  must  be  of  the  font 

*t  -  xi  =  y o  Cf(y0fl  (B3) 

(where  [f(y0)]  /  0,®)  in  order  to  satisfy  the  restrictions 

V  s  0 

Jo 

that  at  the  wing  tip 


9*111 
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A 


(a)  the  chord  is  zero 


(xt  -  x7)  =  0 

t  1  =  0 

(h)  the  slopes  of  the  leading  and  trailing  edges  are  neither 
equal 

AK  ~  xz) 
dyD 

nor  in  the  free -stream  direction 


~d(*t  -  xz) 
dyo  . 

Substituting  equation  (B3)  into  the  denominator  of  equation  (B2) 
yields 


2 

HL  .  _  (B4) 

dyo  v/j0  [f(j0)j 

Equation  (B4)  indicates  that  dr/dyQ  is  singular  of  order  l/2  at 
yc  =  0  for  the  wing  tip  of  figure  6. 
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(x,y»s) 


(a)  Line  vortex  from  (x1>y1). 


(x,y,z) 


(b)  Line  vortex  from  (x2,y2)» 


(c)  Line  vortex  from  (x^,y^)  to  (X2>y2)* 


Figure  2«  -  Superposition  for  obtaining  upwash  induced  by  line- 
vortex  segment  from  x1,y1  to  Xj2»y2« 
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(a)  Infinite  line  vortex  inclined  supersonically 
to  free  stream. 


(t>)  Infinite  line  vortex  inclined  subsonically 
to  free  stream. 


U 


r 

x,xD 


Figure  3.  -  Typical  line  vortices. 
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Figure  4.  -  Lines  of  constant  potential  for  rectangular  and 

triangular  wings. 
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(b)  Elliptic  load  distribution. 

Figure  5.  -  Upwash.  in  z  *  0  plane  an  infinite  distance  behind,  wings 
of  triangular  and  elliptic  loading. 
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Figure  7,  -  Bent-lifting-line  approximation  for  sweptback  wing 
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(a)  Bent-lifting-line  representation  of 
triangular  wing. 


(b)  Unbent-lifting-line  representation 
of  triangular  wing. 


/ 


/F 


/ 


/ 


/ 


/ 


/ 


H) 


/ 


/ 


/ 


\ 


\ 


X 


\ 


\ 


r 


yjVo  / 


/ 


rU 


\  W~o/ 

- ^ — 7 - 

V 
/  \ 

'  \ 


/t' 


/ 


\ 


\  / — Unbent  lifting  line 

W 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


(c)  Unbent-lifting-line  representation  of  rectangular  wing. 


Figure  10.  -  Lifting-line  representations  of  triangular  and  rectangular 

wings. 
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Figure  13.  -  Downwash  along  line  y*0,  z  =  O  behind  triangular  wing  using  unbent  (horseshoe) 

lifting  line  at  3/4 -chord  point.  M  « 
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Figure  14*  -  Downwash  along  line  y»0,  z * 0  behind  rectangular  wing  using  unbent  (horseshoe) 

lifting  line  at  l/2-chord  point.  M  « 


